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Abstract. We have studied the free expansion of a Bose condensate in which 
both the usual s-wave contact interaction and the dipole-dipole interaction 
contribute considerably to the total interaction energy. We calculate corrections 
due to dipolar forces to the expansion of such a condensate after release from a 
trap. In the Thomas-Fermi limit, we find that the modifications are to lowest 
order independent from the total number of atoms. 
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1. Introduction 

In the atomic Bose-Einstein condensates (BECs) realized so far, the atoms interact 
essentially only at very short distance (Van-der Waals interaction) compared to 
the separation between atoms. The only relevant parameter that characterizes the 
interaction is the s-wave scattering length Q). Recently dipole-dipole interactions in 
atomic BECs have attracted increasing interest §, |, §, |[|, 0, |, § 0, ■ These 
interactions would largely enrich the variety of phenomena to be observed due to their 
long-range and anisotropic character. 

Atoms like rubidium and sodium which have been successfully Bose-condensed, 
posses a magnetic dipole moment, but it is very small. We focus here on another atom 
as a possible candidate for the realization of a dipolar condensate, namely chromium. 
Chromium atoms have a larger magnetic moment m = where /is is the Bohr 

magneton, and interact with non negligible dipole dipole forces. These anisotropic 
forces compete with s-wave scattering. If the s-wave scattering length (unknown for 
chromium) is not too large, then dipolar forces could be observed directly through 
small but appreciable modifications in the ground state density distribution or even 
in the occurrence of instability. 

We discuss in this paper the dipolar effects on the free expansion of a dipolar BEC. 
In many cases, trapped condensates are too small to obtain clear images in absorption 
imaging and therefore the condensate is released from the trap to expand in order to 
obtain larger dimension. Because the dipolar effects in chromium are expected to be 
not too large, it is interesting to maximize the dipolar effects in ballistic expansion by 
tuning external parameters like the trap anisotropy. 
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2. Magnetic dipole interaction 

We consider a dipolar gas of neutral atoms with magnetic moment m oriented along 
the magnetic field 

B(t) = B [cos(/3z + simp (cos(fii)x + sin(Oi)y)] . (1) 

The frequency Q, of the radial component of the magnetic field is chosen such that 
ULarmor 3> f2 ^> uJtrap, where Utrap is any of the trap frequencies. This corresponds 
to the assumption that the atoms are not moving during one rotation of the magnetic 
field and the magnetic moments will follow adiabatically the total magnetic field. The 
resulting long range part of the interatomic interaction is [ |Tl| 

Horn 2 ( 3 cos 2 (p — 1 \ / 3 cos 2 6—1 
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where 6 is the angle between the interatomic separation vector r and the z-axes. Here 
m is the atomic magnetic moment and fj,o is the magnetic permeability of the vacuum. 
The angle tp serves as a "knob" to tune the dipolar interaction in real time (<p may 
have a low frequency time dependence). 

3. Thomas-Fermi mean-field analysis 

Our analysis is based on the zero-temperature (T = 0) mean-field approach already 
used in the context of dipolar Bose gases [||, ||, Q]. Such a description can be 
accomplished through the Gross-Pitaevskii equation |l| for the condensate order 
parameter $ (r, t) (normalized to the number of atoms N) 

where M is the atomic mass, to r (co z ) the radial (longitudinal) trap frequency and 
where the mean- field potential V m / is given by 

Wr) = ^|*(r)| 2 + J dr'|*(r')| 2 t/ dd (r-r'), (4) 

where a is taken to be a positive s-wave scattering length. 

A measure of the strength of the dipole-dipole interaction relative to the s-wave 
scattering energy is given by the dimensionless quantity 

_ fi m 2 M 
Yl-KTv-a 

£dd is defined in such a way that when it is larger than one (edd > 1) a homogenous 
dipolar condensate is unstable against collapse in the Thomas-Fermi limit j(| and in 
the absence of an external potential. 

For chromium the s-wave scattering length ao is not known. If we assume 
that acr is equal to the sodium scattering length (a = 2.8 nm) the dipolar strength 
parameter becomes ej^ = 0.29 (m = 6(j,b, where [ib is the Bohr magneton). In 
comparison, for rubidium and sodium the corresponding values of £dd are much 
smaller. In contrast, for heteronuclear molecules, we expect a much larger £dd- 

The starting point for our analysis is the limit of a large number of atoms, namely 
the Thomas-Fermi limit in which the kinetic energy has a very small effect on the 
ground state density distribution. This physically relevant case represents also an 
ideal situation from an experimental point of view because an observable quantity like 
the modification of the condensate shape is independent of the number of atoms. 




Figure 1. Equilibrium condensate anisotropy parameter re = £r/£z versus the 
trap anisotropy parameter A = u z /u> r . ft is calculated in the Thomas-Fermi 
regime with = 0.29. The three curves correspond to three different values of 
tp: the curve for tp = is solid, for tp = 54.7° is dashed and for tp = 90° it is 
dotted. The inset shows Are = (k(<p = 90°) - n(tp = 0°))/re(i ( 3 = 57°) versus the 
trap anisotropy parameter A. 



4. Condensate expansion 

The expansion of the condensate (as well as its equilibrium configuration) can be 
studied using a time-dependent mean-field variation approach similar to Ref. 0, |l5| 
(see also [|ll] ) by assuming a gaussian wave function as an ansatz 

V(x,y,z,t) = }{ e **. , (6) 

V t,xt,yt,z J r,=x,y,z 

where £ v (width) and j3 v are time dependent variational parameters. It is easy to 
derive the following set of equation for the condensate widths 

(f3 v — l^l 2 -)- The various contributions to the total mean-field energy functional H tot 
evaluated on the gaussian wave function are: (a) the kinetic energy in the Thomas- 
Fermi limit i/kin = {NM/4) + iy + £f J (where we neglect here the contribution of 
the zero-point fluctuations {Nh 2 /AM) + 1/^ + 1 /(%))', (b) the potential energy 

in the harmonic trap Hh = (JVm/4) (uj 2 £, 2 + uj 2 ^ 2 + (c) the mean-field energy 

due to s-wave scattering 

H N2H2a 1 • (H) 
mV27r 

(d) the mean-field magnetic dipole-dipole energy, which for a cylindrically symmetric 
condensate is given by 




Figure 2. Condensate anisotropy parameter k after expansion versus the trap 
anisotropy parameter A. n is calculated in the Thomas-Fermi regime with 
£dd = 0.29. The three curves correspond to three different values of ip: the 
curve for tp = is solid, for ip = 54.7° is dashed and for <p = 90° it is dotted. 
The inset shows (k(<p = 90°) — K,((p = o ))/n(<p = 57°) versus the trap anisotropy 
parameter A. After expansion, the observability is improved in the "pancake" 
geometry (large A). 



, 1 + 2k 2 3«; 2 tanh Vl - « 2 

fiK) = — (1 - (10) 

Here, we have introduced the radial width £ r = £ x = £ y . The argument of / is the 
condensate aspect ratio n = £ r /£z- 

The anisotropy of the static dipole-dipole interaction is directly manifest in the 
modification of the shape of the condensate. In the Thomas-Fermi limit the ratio 

n-^H-k (11) 

of the radial (x 2 ) 1 / 2 and longitudinal width (z 2 ) 1 / 2 of the condensate is independent 
of the number of atoms, because the dipolar energy has the same scaling dependence 
as the s-wave scattering energy. In particular k varies as function of the "knob" angle 
ip. The anisotropy parameter k of the condensate can be studied for an equilibrium 
situation, in which the condensate is harmonically trapped, or after a free expansion. 
Generally, the condensate anisotropy k is a function of the trap anisotropy parameter 

A = UOz/uUr ■ (12) 

In the Thomas- Fermi limit k is a solution of the following equation 

A 2 = ,, 2 l-gdd/(«0-e d d/tf£(Ac) 

1 - £dd/(«0 + \£ddH^{n) 

It is straightforward to see the standard n — A result for a trapped condensate if 
dipolar forces are absent. When e^d <C 1 we can evaluate analytically the variation of 
the aspect ratio n and the variation of the "condensate volume" V = ((x 2 ) (y 2 ) (z 2 )) 1 / 2 

6k 3 , 5/ ... 
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Figure 3. Condensate anisotropy parameter k after a free expansion with 
'inverted' dipolar coupling. When the condensate is released, the angle ip is 
switched from ip = to (p = 90°. k is calculated in the Thomas- Fermi 
regime with e = 0.29 as function of the trap anisotropy parameter A. The 
three curves correspond to three different values of ip: the curve for ip = is 
solid, for ip = 54.7° is dashed and for ip = 90° it is dotted. The inset shows 
(n(tp = 90°) — K(ip = 0°))/k(ip = 57°) versus the trap anisotropy parameter A. 



The condensate's aspect ratio k is a more universal parameter than the overall size 
V because it depends only on the trap geometry and not on the number of atoms. 
In particular the variation of k is easy to measure because the value of k in absence 
of dipolar forces is well know (in the Thomas- Fermi limit, it is given by the trap 
anisotropy A; therefore 5 k = k — A). Moreover it is possible to change the tuning angle 
(f and measure the difference n((p = 90°) — n(ip = 0°). In contrast, the variation of the 
condensate volume V is more difficult to measure because it requires the knowledge of 
the number of atoms, a quantity that is usually difficult to determine with sufficient 
accuracy. From Eq. ( |l4| ) it is possible to see that k is not sensitive to the dipolar 
forces for high and in particular for small values of A because the product k§£(k) goes 
to zero in the two limits n — > and n — > oo. The maximum linear response of k, to 
£dcb i- c - 5k ~ .64£ ( jd J is obtained for A ~ 2. 

Figure 1 shows the equilibrium anisotropy parameter k of the condensate versus 
the trap anisotropy parameter A in the Thomas-Fermi regime, calculated for a 
chromium condensate where a s-wave scattering length of a = 2.8 nm is assumed 
(edd = 0.29). The effect of the dipole interaction is to reduce k for ip = 0° or to 
increase k for ip = 90° with respect to the equilibrium value of k in absence of dipolar 
forces (see Figure 1). Note that by varying tp, the total variation of k can be increased. 
Thus, the influence of the dipolar force becomes more visible. 

Now we consider the expansion of the condensate. Like for the static properties, 
the evolution of the condensate anisotropy parameter is independent of the number 
of atoms in the Thomas-Fermi limit. Figure 2 shows the asymptotic values of the 
anisotropy parameter n(t = oo) after expansion. In a finite time experiment these 
values are approximated by the ratios of the velocities £ r /£z- 




(15) 
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Similar to the equilibrium case the effect of a (weak) dipole interaction is to reduce 
k for ip — 0° (or to increase k for if = 90°) with respect to the values n{t = oo) in 
absence of dipolar forces (see Figure 2). In the absence of dipolar forces n(t = oo) is 
always inverted after expansion. 

According to the results shown in Figure 2, a pancake trap will enhance the 
variation of k with if. Assuming a = 2.8nm, the variation of k is 40%. If £dd is in 
reality smaller than our assumption one should think of some scheme to amplify the 
effect. A simple scheme to do that is the following (see Fig. 3): when the condensate 
starts expanding, the sign of the effective dipolar coupling should be changed by 
changing the angle if from 0° to 90°. With this simple technique the observability is 
further increased. 

The particular case of a dipolar condensate initially trapped in an isotropic trap 
(A = 1) deserves some consideration. When the condensate is released its anisotropy, 
i.e. its shape (cigar deformed for if = 0° and pancake deformed for ip = 90°), is 
conserved during the time evolution. Indeed the set of equations (0) with a general 
time-dependent isotropic harmonic potential frequency wo(i) admits a solution of the 
form £ v (t) — £(i)£H with i] = x,y,z (where £° for instance are stationary solution of 
((jj) with LJo(t) — ujq). Similarly the frequency of the monopole mode in an isotropic 
trap is given by VZujq, where loq is the trap frequency, irrespectively of the dipolar 
forces. 



5. Conclusions 



We have theoretically investigate the static properties as well the expansion of a dipolar 
condensate in the Thomas-Fermi regime. We found that the pancake geometry is a 
more suitable trap geometry to observe dipolar effects, since they are enhanced during 
expansion. 

Note added in proof — We are aware of one similar work on the expansion of a 
dipolar condensate Jl6| ]. 
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